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ABSTRACT 



An apparatus was designed and constructed to ineasure the modulation 
transfer function of the atmospheric medium as a function of wavelength 
of light. The modulation transfer function (liTF) in this case was 
considered to be the ratio of the modulation of an optical target at: 
a fixed range to the modulation at a range short enough so that atmo- 
spheric attenuation was negligible. 

Modulation measurements were made using a telescope which focused 
the image of the target on a slit. The image was mechanically scanned 
across the slit by a rotating mirror, and the intensity at the slit was 
measured by a photodetector. 

The preliminary results indicated a linear decrease in modulation 
o 

of 0.1 per lOOOA. 
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I, INTRODUCTION 



With the increasing complexity of optical systems in recent jrears,. 
it has become necessary to find a general figure of merit by which to 
compare them. By viewing the object of interest as an input, and the 
image formed as an output, it has become common practice to use the 
theory of the transfer function as the descriptive link between the 
two. 

Generally speaking, the quality of an image formed by an optical 
system is a function of three things: 1. The total abberation in 

the ray optics design of the system, 2. The (inherent) diffraction 
caused by the wave nature of light, and 3. Manufacturing inaccuracies. 
In the case where manufacturing and abberation effects can be consider- 
ed negligible, the system is said to be diffraction limited. The dif- 
fraction limited case is the closest to ideal that can be reached. 

The observable effect of diffraction is a reduction in contrast which 
increases with the spatial frequency of the target. 

When dealing with the intensity distribution of a periodic target, 
an important comparative parameter is the modulation, defined as the 
difference between the maximum and minimum signal amplitude divided by 
the sum of the two. This parameter is the ratio of (tx^rice) the ampli- 
tude variation to (twice) the average intensity. 

In the most common usage, the plot of modulation versus the spatial 
frequency of the target (cycles per distance on the target) is called 
the modulation transfer function (MTF) curve. Since this is normalized 
to 1.0 at zero cycles, it will be defined in this paper as the ratio 
of the modulation at a fixed distance to the modulation at approximately 
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zero distance for a fixed spatial frequency at a particular Wcivelength 
of received light. Since the atmosphere is the common medium of 
propagation for most optical systems, it is necessary for both, pre- 
diction and evaluation of system performance that an understandings of 
its characteristics be applied to MTF calculations. 

A. THEORY OF MTF FOR A DIFFRACTION LIMITED SYSTEM 

The response of an optical system is determined for basic patfems,. 
such as points, lines, or edges. The responses to these patterns are 
called the point, line and edge spread functions respectively. For 
responses to complex objects, such as periodic targets, the optical 
transfer function is used if both amplitude and phase are of concern 
and the MTF if amplitude only is desired. Since any object may be 
considered a collection of points, there is a connection between the 
point spread function and the transfer functions. 

The point spread function is* the normalized diffraction pattern 
formed by a point object. It is usually plotted as illuminance as a 
function of distance in the image plane. 

Once the point spread functions for a sensor system are known, the 
illuminance of any object can be obtained by summing the point spread 
functions of the target object. This process can be represented by 
a convolution integral. As a general example, consider an object plane 
(x^) and an image plane (x^). Let the function spread 

function SCx^) define an object. The illuminance of the image, I(x^), 
,is given by the convolution of the object functions: 

oo 

I(Xi) = ^ s(.x^)0(x^-x^)dx^ 

~ OO 
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Figure 1 shows the convolution of an object function and a spread 
function to produce an image function. The image function does not 
exactly reproduce the object because only three points in the image 
function are illustrated. The appearance of the argument (x^-x^) 
implies that one function is translated or scanned along the other 
function so that each point in the object is traced through the 
system to the image, and the spread function is then scanned over 
the image point. 

For a linear system, the response due to a number of inputs is 
the product of the transfer functions of each. In an optical system, 
the mathematical foundation is the Fourier transform: 

oo 

G(w) = J"g(x)exp(-iwx)dx 

— oo 

In any real system, the spatial variable is usually at least 
dimensional. Taking the Fourier transform of both sides of the con- 
volution integral , 
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iCw) = 



oo 

S(x2)0(Xj^-X2)exp(-iwx^)dx2dXj^ 

-oo ®o oo 

( SCx2)dx2 0(x^-X2)expC-£wx^)dx^ 

^ —00 

— oo 

r /” 

I S(x^)dx^ j 0(u)exp(-iwxL)expC-’iwx2)du 



where u = 

L 2 



= ^S(x2)exp(-iwx2)dx2 y^O(u)expC.-iwa)da 



So that I(w) = S(v:)0(w). 



The spread function can then be written S(w) = . When, normalized, 

the Fourier transform of the spread function is the optical transfer 
function. 

S(x)exp(-iwx)dx 

La 

OTF(w) = 



oo 

/ 



•CO 

S(x)dx 
w 



The MTF as a function of k = . can be obtained by: 

27T 



MTF(k) = 



OTF(^) 



It should be pointed out here that k is the spatiail frequency in lines 
\ 

per distance, and that w is therefore a spatial angular frequency, in 
radians per distance. Recalling that 

exp(-iwx) = exp(-i27rkx) = cos( 2Trkx)-isin( 2Trkx) , 



MTF(k) = 



00 00 
^ S(x)cos(2TTkx)dx-i ^ S(x)sin( 2irkx)dx 

I j S(x)dx 



The numerator is the vectorial sum of the cosine and sine Fourier 
transforms of S(x). 
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As an example of the calculation of the MTF of a spread. func*tion 



consider a function defined as: 



SCx) 



1. ^<x<| 

0 5 elsewhere 



This might represent the spread functon of a scanning spot of diameter 

a with spatial frequency — . Since the spread function is symmetrical,. 

a 



I 



a/2 

S(x)sin(2Trkx)dx = 0 for all, k 

a/2 



-a/2 

a/ 2 f 

I S(x)cos(2TTkx)dx = 2 I cos(2Trkx)dx 



-a/2 



-a/2 



= 2 



sin( 2Trkx) 



a/2 



^ 2Trk 

sinCirka) 

Trk 



-1 



a/2 



The denominator, j S(x)dx = a, so that 
-aV2 

MTFCk) = 



sin(Trka) 



irk a 



It can be seen that the MTF vanishes at integer multiples of ~ . 

The existence of the MTF beyond its cutoff frequency of — is 

a 

referred to as spurious or false resolution. This can be explained 
by considering a bar chart of lines and spaces of widtli a. If the 
scanning spot of diameter a is positioned so that the center of the 
spot traces the edge of each successive line, the output of each 
scan v;ill be the same because each will scan the same amount of line 
and space. If the spatial frequency is increased, the lines and 
spaces become narrower, and therefore closer together. In this 
case, the amount of line and space covered each scan varies from 
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nearly all line to the opposite, and the varying oritput looks, like 
resolution, although no useful information can be obtained from LL* 
In the case of a diffraction limited lens, the spread function 
is the pattern consisting of the central Airy disk surrounded by 
progressively fainter rings. The MTF is the Fourier transform of 
this diffraction pattern. Consider the general case of Fraunhofer 
diffraction of a circular aperture. Figure 2 illustrates the co- 
ordinates used. 

/ 

e 

2 



Figure 2. Fraunhofer diffraction of a circular aperture. 

It is desired to find the illumination at a point due to a 
source at point P after the radiation has passed through an. aperture 
A. The amplitude of P^ is found by summing all the contributions 
from the area elements in A. 

U(P' ) = C r U(X.Y) dA = C C f °>=P(iHrtr-))dA 

Because kr is large, the numerator changes rapidly with r, while the 

denominator changes slowly over the range of integration, so long 

as rr’ >> X , Y 

max max 
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« 






■r 



r/M 






= CCx-X)^ + Cy-Y)2 + z2] 



r/2 



= r 



(1 + -2xX - 2yY^ 



1/2 



:: r - 
o 



xX + yY + + Y^ - CxX + yY) ^ 



2r 



2r 



A similar expression can be obtained for r^. In these expressions,, r^ 
and r^ are the distances from the origin to P and P* respectively. Now 
let 1 and m be the first t>ro direction cosines for P and 1* and m*’ be 
those for P' . These cosines may be written: 



- “X -y - , X* V*' 

1 ~ , m - — 5 I “ — r * m^ = 



Substituting, 



)Y 

Z 



U(P’) = Cexp(ik(r^ + f f fik((l - l' )X + (m-ffl'- 

"a ’ 

if (1, t 1^) (Xl Yl) - - (1.'.^.^.“'..^)"] 

2/^rr*^ r r’ d 

(X^ + Y^) 



dXdY 



In the- case where 



max, the terms in the second 



bracket may be neglected. This is the Fraunhofer diffraction approxi- 
mation. The amplitude expression may be simplified to: 

U(Al,Am) = C exp[ik(XAl + YAm)]dXdY 
A 

Assuming the aperture does not transmit perfectly,, and calling its 
transmittance T(X,Y), 

00 

U(Al,Am) = C jy T(X,Y)exp(ik(XAl + YAm))dXdY 
The amplitude of the diffracted vrave is the t>70-dimensional Fourier 
transform of the aperture transmittance, provided the aperture co- 
ordinates are measured in units of ~ = . Using the Fraunhofer 

K 2tt 
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approximation, the linear coordinates may be writteriy x = r^Al,. 



y - r*Am 
^ o 



U( 



x,y) = c fj 



TCX,Y)expC^ (,xX -r yy))dXdY 

00 O 



Consider a circular aperture of radius a. In order to integrate 
in polar coordinates, we substitute 

X = Rcos0, Y = Rsin0 

X = rcos<{), y = rsin<() 

dA = RdRd0, f = r' 

’ o 



a 27 t 



U = c ^ exp J ■ - (cos0cos(f) t sin0sin<j))J 



d0dR 



o o 
a 2tt 



U = c 



^ ^ j exp cos(0-4)) d0dR 



^2tt 

Recall that I exp(ixcosa)da = 2ttJ (x), where J (x) is the zero order 

Bessel function. Since the integrand is periodic, we may add any 

finite amount to both limits simultaneously without changing the 

/ a 

RJ ( — ^)dR. The recurrence relation for 
o r 

Bessel functions states that 



-.(x)) = (x) , so that 

dx n+1 n 



/ 



xJ (x)dx = xJt(x) 
o 1 



U = (ira^c) 



kra' 



The irradiance is the square of the amplitude, or 
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I 



kra 



I • where I 
o o 



AP 



A = aperture area, P = total power 

The corresponding distribution is called an Airy patlem, and the 
area enclosed by the first miniinum is an Airy disk. The diameter of 
the disk turns out to be approximately 1,22 — ^ , where is the 
apex half angle of the cone converging on the image points 

To determine the amplitude spread function, consider Figure 3- 
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The spherical wavefront of radius R has its vertex. A, af the 
center of an aperture, Q, which can be considered the exit pupil of 
a diffraction limited lens. P(x,y) is a point in the image plane.. 

The amplitude of the wave at P is obtained by integrating from all the 
area elements dx^dy^ over the wavefront. For each point B'(x’ ,y') in 
the aperture, a projection is made to a point B(x\y') on the wave- 
front. A wavefront area ds corresponds to each dx^dy^* The pro j" action 

ds “5 

is made such that - 5 — — r = cos'^G. The amplitude U at P cart be 

dx My ' ^ o 

obtained by integrating over the area elements at Q. 

rr ^ 

Uo(x,y) -j\ [Eexp(iks)cos^0 dxMy’] 

$ r/ -ik 

~ exp(ikR) H expC(— + yy'')] dx'dy' 

E = aperture illumination 
^ = total energy flux in the aperture 
Sq = aperture area 
XV 

Assumptions: "w cosO = 1; uniform illumination* 

R R 



The integral above is the Fourier transform of the anerture function. 

2n 

For a circular aperture this becomes U (r) = ; — , where k - r — , 

^ o krA A 

o 

X is the wavelength in vacuum. A is. a numerical representation of the 

1/2 

aperture, and r = (x^ + y^) . Note that this is the square root of 

the Airy function. A pupil function can be defined, V(x\yM = 

T(x’ ,y * )exp( ikv(x’ ,y ’ ) ) , where T is the amplitude and v the phase 
deviation in the aperture. VJriting the amplitude spread function in 
terms of the pupil function. 



U(x,y) 
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The amplitude spread function is, then, the inverse Faurier 
transform of the pupil function. 

To obtain the transfer function of a lens, it is necessary to 
take the Fourier transform of the spread function, llew variables 

X ^ V ^ 

are defined; v' = — , and v' = ^ , so that 
X. AK y AK 



oo 

U(x,y) = C )/ VU V ’ ,v’ )expt”i(xv’ + yvOjdv dv 

JJ x y^ x^y"^x3 






V(x’ ,y ’ )dx’dy ’ 

This is the inverse Fourier transform of V’, the pupil function,, 
written in the new variables. For an incoherent object, the transfer 
function is the Fourier transform of the intensity spread function,. 



rp/ N ^ [E(x,y)] ^ 

g ;CE(x,y)] ■ = 



U(x,y) 



E(x,y) = C \V)] Gh \v’)] 



Applying the convolution theorem of the Fourier transform, 

c^{Eh-kv)] 



T(v ,v ) = 
X y 






, but 



S[{ifi'*[v(x):)" = v''(-x) 

Thus the convolution changes to correlation, 

/ oo 

f(x’ )g(x’-x)dx’ .. 

u 

The denominator ,^^(E) = (0,0) = CT, the total flux through 

the aperture. V/hen normalized, T(v ,v ) = , vxhich represents 

x’ y ^ " 
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the value of integrated over the common area when two figures 

of identical cross section are displaced by Taking the example 



of Figure 4, where tvro circular apertures are shown, the ITTF 




CGinmon area 
total area* 



Figure 4. Correlation of two identical circular apertures.. 

V V V 

• • A X XX 

From Figure 4 it can be seen that rcos9 = -r— , so that cos<J> = -r - — = — — . 

2 2r D 

2 

The common area is given by (rsin(J))(rcos(})) J = 2r^( (})-sin(|)cos<f5-) 



MTF = — ( (j)-sin<{>cos(})) 



If the maximum spatial frequency is k^, then the normalized spatial 



frequency = — = — = cos(}), and sin(j) is equal to (l-cos^(|)) 



1/2 



MTF 



2 f -Ik . k ,k 



This is the Fourier transform of 



2J^(^kQX 



irk X 
o 



1 



the Airy pattern, so 



that the MTF is that of a diffraction limited lens. The interesting 
point is that the MTF is the convolution of the aperture with itself. 



B. THEORY OF THE EXPERIMENT 

The wavelength dependence of the MTF of the atmosphere was to be 
obtained by optically scanning a 4' x 6* bar chart made of alternating 
black and v:hite strips. Each strip was ten centimeters wide, with a 
resulting spatial frequency of five lines per meter. The angular scan 
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■t 






provided by a rotating mirror was a total of 0.003 8 radians y and was 
interrupted by a mechanical chopper at the rate of 29 chops per second 
to provide a reference zero. The chopper and rotating mirror were 
constructed as shown in Figure 5. The wavelength of the received light 
was to be varied through the visible and into .the infraredy and the 
modulation observed and referenced to the modulation obtained at. 
essentially zero distance to yield the MTF. 

C. EVOLUTION OF THE EXPERIMENTAL APPAPJ^TUS 

The block diagram of Figure 6 was a common factor throughout the 
research effort , but several modifications were made to the optics in 
an effort to improve signal amplitude and accessibility of the slit for 
visual target monitoring. 

The initial optical arrangement was as shown in Figure 7y using 
the monochromator section of a Perkin-Elmer Model 13 infrared spectro- 
photometer. In this arrangement, the light passed through the entrance 
slit, was dispersed by the prism, passed through an array of mirrors 
and slits, and was finally reflected up through an aperture in the 
cover of the monochromator to a silicon photocell by mirror M6. The 
detector circuit was as shown schematically in Figure 8y except that 
the rectifier diode had not yet been added. Wavelength scanning was 
achieved by changing the angle of M2 to cover the entire prism 
dispersion pattern. 

A problem was discovered with this arrangement in that the losses 
suffered by the light in passing through 9 mirrors (including two in 
the telescope), through the prism twice, and into the detector were 
so great that the signal amplitude out of the detector was too low to 
be observed above the noise on the Tektronix 451 oscilloscope. The 
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Chopper and Scanner Mechanisms. 
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Figure 7. Initial Optical Arrangement. 
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f?igure 6. Detector Circuit Schematic 



first effort to improve this situation was to shift from a. Type E 
Low Level Differential AC Preamp plug-in scope module (which had a 
lower range of 50 microvolts per cm) to a type 1A7 High Gain Ilif— 
ferential Amplifier Cwith a lower range of IQ microvolts per cm).. 

This was still unsuccessful. With the entrance slit full, open, at 
2 mm, the white bars on the target were visible on the scope:,, hut 
the black bars were lost in the noise. 

The next modification v;as to bypass the monochromator entirely,, 
and place the detector directly behind the entrance slit.. This gave 
detectable signal amplitude, but a problem of overshoot on the scope 
coupled with capacitive discharge through the detector circuit 
produced a negative response of approximately the same amplitude as 
the signal. At this point, it was decided that at least part of' the 
problem could be solved by dealing with the capacitor. The value of 
the capacitor was changed ttfice, ending with ^047 microfarads . ■ It was 
also observed at this time that the photocell output as observed on 
the scope v;as overshooting in the negative sense on. each chop, either 
as a result of the internal oscilloscope or detector circuit charac- 
teristics. 

By partial dismantling of the system, it was possible to determine 
visually that the chop on each cycle was occurring while the image 
of the target had been sv;ept past the slit, which implied that there 
was only background present when the chop occurred.. This problem was 
corrected by rotating the cam on its shaft. 

At this point in the research , an attempt vras made to change 
detectors from the silicon photocell to an InSb photovoltaic, liquid 
nitrogen cooled detector manufactured by Barnes Engineering. In orde: 



22 



to accommodate the detector, shown in Figure 9, the optics were 
rearranged as shown in Figure 10. The documents pertaining to 
the detector failed to shovr the associated circuitry required,, and,, 
as a result, no useful data were obtained from it- An attempt was 
made to take the output directly from the detector to the scope,- 
since it operated at zero volt bias and had low maximum current and 
voltage ratings, but no signal was observed. 

In order to alleviate the capacitor discharge problem, the 
rectifier and resistor were added to the detector circuit so that 
any discharge would occur through the internal circuit vice through 
the output. 

Since phone calls and local efforts x^rere unsuccessful in getting 
results from the InSb detector, the silicon photocell was reinstalled 
using the optics as shovm in Figure 10. Wavelength variation was 
achieved through the use of interference filters. The results obtained 
will be presented in Section II of this paper. 
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Figure 9. InSb Detector 
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;ure 10, Path of Lif^ht thi*ough Telescope and Optics 



II. RESULTS OF THE * EXPERIMENT 



The following data were obtained using a target illuTnina1:ion of 
approximately 200 watts per square meter by incandescent light, a 
chopping frequency of 29 per second, and essentially zero attenuation 
due to distance. 

The reference modulation of the target based on 46 measurements 
at distances between 20 and 120 meters was Q..647 + .011.. This sirould- 
be normalized to 1.0 in dealing with future data. Graph 1 is a plot, 
of modulation versus wavelength at zero range. The reliability of 
this graph is questionable since the interference filters used were 
of varying transmittance. In order to resolve this problem,, a sel 
of matched, high transmittance filters should be used. It is expect- 
ed that the trend of this graph will reverse to some extent at the. 
higher wavelengths, using an infrared detector at greater range, due 
to the infrared window of the atmosphere. 

Figure 11 shows scope • photographs of typical signal presentations.. 



26 



MODULATION VS WAVELENGTH AT ZERO RANGE 



/ 



o 



/ 



c 

c 

c 

no 



/ 



© 

/ 



I 



c 

o 



/ 

I 



c 

c 

c 

r- 



/ 



/ 



I 



o 

o 

LfN 

VO 



/ 



/ 



o 

c 

o 

VO 



I ® 



© 



© 



/ 



/ 

o 



o 



o 

/ 

t I ; 

O m 

LT\ C\i 

• • « 



o 

m 



c 
- c 
o 

LTs 



o 

o 

LTV 



o 



Noiivinoo/i gszit:v7ihom 



27 



WAVELENGTH OF RECEIVED LIGHT IN ANGSTROMS 




Figure 11a. Typical 
oscilloscope presentation. 
Highest amplitude is white 
bar, double peaks are black, 
and chop is between. 




Figure lib. Expanded 
scope trace. 
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